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ABSTRACT 

The Regional Acid Deposition Model (RADM) 
hns been applied to several of the field experiments which 
were part of Ihc Acid Models Operational and Diagnostic 
Evaluation Study (Acid MODES). The experiment which 
was of particular interest with regards to ozone photo- 
chemistry involved horizontal zig-zag flight patterns 
(ZIPPER) over an area from the eastern Ohio River valley 
to the Adirondacks of New York. Model simulations by 
both the standard resolution RADM (Ax = 80 km) and the 
nested grid RADM (Ax = 26.7 km) compare well to mea- 
surements in the low emission regions in central 
Pennsylvania and upstate New York, but underestimate in 
the high emission upper Ohio River valley. The nested 
simulation docs considerably belter, however, than the 
coarse grid simulation in terms of horizontal pattern and 
concentration magnitudes. Analysis of NO x and HO x 
concentrations and photochemical production rates of 
ozone show that the model’s response to large point source 
emissions is very unsystematic both spatially and tempo- 
rally. This is due to the models inability to realistically 
simulate the small scale (subgrid) gradients in precursor 
concentrations in and around large point source plumes. 

1 . INTRODUCTION 

Photochemical Eulcrian grid models arc being in- 
creasingly relied upon for understanding, assessment, and 
regulation of mcso-scalc and regional ozone problems. 
Therefore, it is extremely important to understand the real- 
ism and limitations of these models. It is often assumed 
that if grid resolution is made sufficiently small the model 
can realistically simulate the complex interactions of pho- 
tochemistry and dynamic transport and dispersion pro- 
cesses. This is true only where the grid can sufficiently 
resolve the gradients of reacting chemical species. As 
modeling emphasis shifts from small urban areas to larger 
regions (c.g. Southern Oxidant Study, Lake Michigan 
Ozone Study, San Joaquin Air Quality Study), models with 
coarser grid resolution arc being applied to photochemical 
problems. These modeling efforts arc complicated by the 
presence of very large point source plumes which dominate 
NO x emissions in many areas. Ozone formation associ- 
ated with power plant plumes has been observed many 
times (Davis cl al., 1974; Miller ct al., 1978; Gillani and 
Wilson, 1980). Since ozone photochemistry is a notori- 
ously non-linear process, the ability of Eulcrian grid 
models to simulate ozone production in the vicinity of 
large point source plumes, where gradients of precursor 
concentrations arc particularly severe, needs to be as- 
sessed. 

2. DESCRIPTION OF AIRCRAFT MEASUREMENTS 

The specific measurements involved in this study 
were made by two specially instrumented aircraft on August 
31, 1988 as part of AcidMODES. Two types of flight 
patterns, designed to diagnose various aspects of RADM 


performance under clear sky conditions, were flown in a 
corridor extending from the eastern Ohio River Valley to 
northeastern New York. This study focuses on the ZIPPER 
flights, a horizontal zig-zag pattern flown within the 
mixed layer (about 1300 m AGL), in order to discern hori- 
zontal spatial patterns and gradients. Both aircraft 
recorded continuous measurements of gaseous SO2, NO, 
NO2, NOy, O3, H202> and CO as well as ambient tempera- 
ture, moisture, pressure, and winds. 

3 DESCRIPTION OF MODEL 

The Regional Acid Deposition Model (RADM) is 
a comprehensive Eulcrian atmospheric chemistry grid 
model which includes statc-of-scicnce representations of 
gas-phase chemistry, advective transport, subgrid vertical 
mixing, dry deposition and cloud processes including 
aqueous chemistry, convective mixing, and wet deposition 
(Chang ct al., 1991). The modeling system incorporates 
hourly emission inventories, both anthropogenic and 
biogenic, land-use information, and gridded hourly meteo- 
rological fields provided by simulations of the Mesoscalc 
Meteorological Model, Version 4 (MM4) with 4-D data 
assimilation (Stauffer ct al., 1990). In addition to the 
standard regional scale version of RADM which has a hori- 
zontal grid cell size of 80 km, a nested grid version of 
RADM (Plcim et al, 1991) which increases horizontal 
resolution by a factor of 3 (Ax = 80 km/3 = 26.67 km) is 
used in this study in an attempt to investigate the effects of 
model resolution on mesoscalc ozone simulation. A revi- 
sion to the model since the description by Chang ct al. 
(1991) is the addition of a non-local closure scheme for 
subgrid vertical transport in the convective boundary layer 
called the Asymmetric Convective Model (Plcim and 
Chang, 1992). 

4. COMPARISON OF MEASURED AND MODELED OZONE 
CONCENTRATIONS ON A REGIONAL SCALE 

Figure 1 shows the horizontal spatial distribu- 
tion of ozone measured along the ZIPPER flight path (zig- 
zag line in Fig 1) and interpolated between the flight legs 
using a simple inverse distance weighting technique. In 
order to facilitate comparison with model simulations, the 
measured and interpolated values arc spatially averaged 
onto the same grid used by the nested grid model (Ax = 
26,7 km). Direct comparison to model simulations is 
complicated by the fact that the ZIPPER aircraft measure- 
ments took almost 4 hours to complete (1530 - 1920 
GMT). Model simulations of this period, however, showed 
only small changes with time. Therefore, ozone concen- 
trations at 18 GMT as simulated by the nested and coarse 
grid models arc shown in Figure 2 for comparison to Figure 
1. Model level 5 (~600-900 m) was used for comparison to 
the aircraft measurements even though the aircraft’s alti- 
tude was mostly in level 6 because the current model gener- 
ally underestimated the height of the convective mixed 
layer and level 5 is the highest model layer which was 
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entirely within lire model’s mixed layer for the duration of 
the flight. 



Figure 1. Airborne ozone concentration (ppb) measure- 
ments along ZIPPER flight path and interpolated between 
legs. Data is averaged onto a grid with 26.7 km grid cells 
(same as nested grid model). 



Nested Grid 



Coarse Grid 

Figure 2. Ozone concentration (ppb) simulated by the 
nested RADM (Ax = 26.6 km) and coarse grid RADM (Ax = 
HO km) in level 5 (~600-900m AGL) at 18 GMT. 


Both nested and coarse grid model results (Figure 
2) show very good similarity in terms of both spatial pat- 
tern and concentration magnitude to the measurements in 
the low emission regions in tire northeastern portion of 
the ZIPPER. Throughout central and eastern PA and NY 
both model simulations arc mostly within about 5 ppb of 
the measurements. In the higher NO x emission region 
near the eastern Ohio River valley, however, the two simu- 
lations diverge significantly from each other and the 
measurements. For example, the location of the pro- 
nounced depression in the ozone concentration Field near 
the OH-PA-WV border was simulated with much better 
accuracy by the nested model but the depth of the hole was 


more closely simulated by the coarse grid model, although 
some 80-90 km to the south. The ozone maximum mea- 
sured at the southmost ZIPPER comer near the Ohio river 
and into West Virginia (Figure 1) shows a broad area of 
concentrations in the 80-85 ppb range with a peak at about 
85 ppb. The nested model simulation shows peak concen- 
trations in this same area of 70-75 ppb at 18 GMT while 
the coarse grid simulation shows concentrations of 55-60 
ppb in this area increasing to 60-65 ppb to the southwest 
(downwind). In general the nested model shows a marked 
improvement over the coarse grid simulation in the high 
emission Ohio River Valley region but is still unable to 
simulate lire peak measured values. 

5. ANALYSIS OF OZONE FORMATION IN THE OHIO 
VALLEY REGION 

In this section the high ozone concentration 
region encountered by the aircraft near the southern end of 
the last comer in the ZIPPER flight path is further exam- 
ined in an attempt to understand how the model simulates 
ozone formation in the vicinity of very large NO x emit- 
ting elevated point sources. This section of the flight 
(~100 km long), where the aircraft crosses the Ohio River 
into WV and then doubles back into Ohio, is shown as a 
time series of 5 second average ozone and NOy concentra- 
tion data in Figure 3. The two peaks in the NO y 



Figure 3. Time series of aircraft measurements of ozone 
and NOy during ZIPPER flight when the G-l aircraft flew 
through the same plume in western WV twice. The vertical 
line through the middle of the plot indicates where the 
aircraft made a 180° turn. 

concentration which coincide with sharp holes in the 
ozone field arc clear indications that the aircraft flew 
through freshly emitted plumes. In fact, since the aircraft 
turned around at the southern comer of the ZIPPER (western 
WV) in the middle of this time series (indicated by the 
vertical line in Figure 3), the two NOy peaks are actually 
the same plume sampled twice. Using data from the navi- 
gational system onboard the aircraft the source of the 
plume was identified as a power plant which emits about 
50,000 tons/year of NO x located in southwestern West 
Virginia near Charleston. During the first pass the aircraft 
was about 3 km from the source and about 6 km on the 
second pass. The location of the ozone depression corre- 
sponds very closely with the NOy peak due to the dcplc- 


71 





lion of ozone by reaction with freshly emitted NO in the 
plume. The peak in ozone concentration, up to 105 ppb 
on the first pass, immediately to the northwest of the 
ozone hole indicates photochemical production of ozone 
at the periphery of the plume where ambient air with VOCs 
and radicals is rapidly mixed with the NO x emissions. The 
smaller peak in the NOy profile which coincides with the 
ozone peak in the first pass through die plume suggests 
that tins is a more dilute portion of the plume which broke 
away from the main plume and mixed more quickly with 
background air. 

These measurements arc an example of photo- 
chemical ozone production occurring very close to die 
source in large point source plumes. Plume photochem- 
istry can be understood in terms of three stages (Gillani 
and Wilson 1980). The first stage is characterized by very 
concentrated NO causing a distinct depression in the back- 
ground ozone field. In stage two, ozone starts to be pho- 
tochcmically produced at the plume edges where the NO x 
has been sufficiently diluted with background VOCs and 
radicals to not only cause ozone levels to recover their 
background concentrations but to exceed them. A plume is 
in stage three when die ozone depression in the plume 
center, which still exists at stage two, is no longer evi- 
dent. The aircraft measurements shown in Figure 3 seem to 
be an example of the early phases of a stage 2 plume since 
die plume core still shows a significant ozone deficit while 
die edges, particularly the north side of the first pass, 
show ozone excess. 1 

The simulation of plume photochemistry by 
Eulcrian grid models is quite different from the three stages 
described alx>vc. The evolution of plume chemistry is 
closely coupled to plume dispersion, however, this entire 
process is subgrid to models which have grid cell sizes on 
the order of tens of kilometers. Therefore, by examining 
the nested model simulation in this same area we hope to 
learn more about the model's response to plume photo- 
chemistry. Figure 4 shows several chemical quantities as 
simulated by the nested grid model at 17 GMT. The first 
plot (top) shows NO x concentration at level 5 (-900 m). 
NO x is clearly dominated by urban areas and the many 
coal -fired power plants along the upper Ohio River. The 
second plot (middle) shows HO x (H0 t H0 2 ) concentrations 
vc Inch very closely anti-correlate with NO x concentration. 
This is because the radical termination reactions are pro- 
portional to NO x concentration while radical initiation 
reactions arc not, thereby resulting in net radical depletion 
where NO x concentrations arc high. For this discussion, 
the HO x field is used as a surrogate for total radicals. 
Similarly, the rate of reaction between NO and H0 2 , which 
is shown in the third plot (bottom), should be a good 
surrogate for photochemical ozone production by both 
H() 2 and R0 2 reactions. This plot shows peaks of ozone 
formation rales ir. regions of moderate NO x concentration 
such as southwestern WV and near Washington DC. Where 
NOx is highest, such as near the OH-PA-WV border and the 
western tip of Lake Erie, the radicals have been terminated 
to such a degree that the NO-radical reactions arc relatively 
slow. Therefore, ozone formation is slower than ozone 
titration by NO emissions resulting in relative minima in 
Lhc ozone field. Near the plume in southwest WV, on the 
other hand, NO x concentrations arc considerably lower and 
HO x concentration much higher resulting in a favorable 
mix for ozone production. Note that, die greater amount of 
ozone formation in the Washington area (see the ozone 
plots in Figure 2) may be due to the greater amounts of 


VOCs in the urban plumes and the more diffuse nature of 
the NO x emissions (mostly area rather than point sources). 




Figure 4. Nested model simulations in level 5 at 17 
on August 31, 1988 of NO x concentration (lop), 
concentration (middle), and rate of reaction between 
and NO (bottom). 


In order to further investigate the models behav- 
ior, a set of time series for the grid cell in southwest WV 
where the nested simulation shows an isolated peak in 
ozone formation arc shown in Figure 5. The plots arc for 
model level 1 (-0-80 m) where hours 0-31 correspond to 
00 GMT 8/31/88 - 07 GMT 9/1/88. The daytime rise in 
ozone concentration of about 40 ppb from hour 1 1 to hour 
13 is coincident with a peak in the reaction rate between 
NO and H0 2 at 12 GMT of about 20 ppb/hr thereby show- 
ing the ozone increase to be primarily photochemical. 
The reason for this sharp peak in ozone formation can be 
understood through examination of NO and HO x concentra- 
tions. NO concentrations begin increasing at 1 1 GMT (7 
am EDT), probably from the morning rush hour in nearby 
cities such as Charleston, WV, while HO x concentration 
initially rises in the morning, peaking at 12 GMT, then 
plunges to a minimum at 15 GMT. Therefore, the rapid 
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ozone formation at 12 GMT results from the combination 
of the onset of significant photochemistry, causing an 
increase in HO x concentration, and increasing NO emis- 
sions. These favorable conditions, however, arc short 
lived since the continuing rise in NO concentration leads 
to a decrease in HO x , and by analogy other radical species. 
Thus, the phenomena of isolated peaks of ozone produc- 
tion, which was seen in the spatial model fields, is also 
evident in the temporal dimension such that ozone produc- 
tion during a single hour can dominate when conditions are 
most favorable. 

6. DISCUSION AND SUMMARY 

In general, whether a particular grid cell contain- 
ing a large point source plume results in an ovcrcstimation 
or underestimation of ozone production depends on a vari- 
ety of factors such as emission rate, grid cell size, disper- 
sion (wind speed and mixing height), and chemical back- 
ground (VOCs and radicals). Thus, for the same condi- 
tions, models of different grid resolution may result in 
quite different simulations. The problem is the grid mod- 
el’s inability to simulate the subgrid processes of plume 
chemistry. The particular ease discussed here suggests that 
increasing grid resolution from 80 km to 26.7 km im- 
proves the simulation of ozone formation. The highest 
ozone concentrations observed by the ZIPPER flights in 
the area downwind of the large point source plumes in die 
eastern Ohio Valley, however, were significantly underes- 
timated even by the nested model. Analysis of the nested 
model simulation of this area shows that just a few grid 
cells, which happened to have a favorable mix of precur- 
sors, were responsible for the bulk of the ozone formed in 
this region. In addition, the time series in Figure 5 show 
that in these grid cells ozone formation occurred very 
sporadically when changing concentrations were briefly 
favorable for ozone formation. To improve model perfor- 
mance subgrid plume chemistry should be simulated either 
by a reactive plumc/puff model or through parameteriza- 
tion of subgrid chemistry based on emissions, chemical 
and meteorological conditions. 
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Figure 5. Time series (00 GMT 8/31/88 - 07 GMT 9/1/88) 
of nested model simulations in the grid cell corresponding 
to the plume encountered by the G-l aircraft shown in 
Figure 3. From the top, the plots show 03, NO, rate of the 
reaction between HO2 and NO, and HO x in model level 1 
(-0-80 km). 
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